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ABSTRACT 

During its 2003 outburst, the black-hole X-ray transient H1743-322 produced 
two-sided radio and X-ray jets. Applying a simple and symmetric kinematic 
model to the trajectories of these jets, we determine the source distance, 8.5 ±0.8 
kpc, and the inclination angle of the jets, 75° ±3°. Using these values, we estimate 
the spin of the black hole by fitting its RXTE spectra, obtained during the 2003 
outburst, to a standard relativistic accretion- disk model. For its spin, we find 
a* = 0.2 ±0.3 (68% hmits); -0.3 < a* < 0.7 at 90% confidence. We rule strongly 
against an extreme value of spin: a* < 0.92 at 99.7% confidence. H1743-322 is 
the third known microquasar (after A0620-00 and XTE J1550-564) that displays 
large-scale ballistic jets and has a moderate value of spin. Our result, which 
depends on an empirical distribution of black hole masses, takes into account all 
known sources of measurement error. 

Subject headings: black hole physics — stars: individual (HI 743-322) — X-rays: 
binaries 



Introduction 



About 50 stellar-mass black holes have been discovered and about two dozen of these 



have been well studied at optical or radio wavelengths (IRemillard &: McClintockll2006uOzel et al, 
20101 ). They are all accretion-powered X-ray sources located in X-ray binary systems. In each 
system, the X-ray source is fueled by gas that feeds from a mass-donor star into the black 
hole's accretion disk. Within a few hundred kilometers of the black hole, the gas reaches a 
temperature of ~ 10^ K and produces a luminosity that can approach the Eddington limit 
(~ 10^^ erg s^^). More than 80% of such sources are transient, with outbursts lasting a 
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year or so followed by years or decades of quiescence. Typically, the host binaries have short 
orbital periods (P ~ 1 day) and are comprised of a low-mass (< 1 Mq) secondary star and a 
~ IOMq black hole. Presently, neither the masses nor the orbital period of our featured sys- 
tem, H1743-322, are known. Nevertheless, as we now describe, the wealth of data available 
for H1743-322 (hereafter H1743) strongly indicates that it is a typical short-period black 
hole transient. 

Studies o f X-ray spectral and timing data leave little doubt that HI 743 contains a black 
hole primary (IKalemci et al.ll2006l : iMcClintock et al.ll2009l : iMotta et al.ll201Cll ). notwithstand- 
ing the lack of dynamical evidence. While large outbursts of H1743 occurred in 1977, 2003 
and 2008, our focus here is on the major 2003 outburst. During this 9-month active period, 
HI 743 was observed 170 times using the PCA and HE XTE detectors aboard the Rossi X- 
ray Timing Explorer (RXTE) (IMcClintock et al.l 120091 ). The source exhibited two distinct 
phases of evolution: (1) During the first three months (when the source was observed on an 
almost daily basis) HI 743 flared continuall y and violently and was in the steep power-law 
(SPL) or intermediate (SPL:Hard) state (see lRemillard fc McClintockl for discussion of these 
X-ray states). On the 47th day of outburst (MJD 52766), an event of central importance 
occurred - the radio/X-ray jets we model were launched during an intense power-law flare 
(discussed below). (2) During the next four months, the source was locked in the thermal 
dominant (TD) state, and the source intensity decayed smoothly and monotonically. It is 
primarily these TD-state data that we use to determine the spin of HI 743. 

An important X-ray timing result derived from the 2003 outburs t was the discovery 



of a pair of quasi-per iodic oscillations (QPOs) at 240 Hz and 165 Hz (IHoman et al.l 12005 



Remillard et al.l 120061 ). Similar high-frequency (HP) QPOs with a commensurate frequency 
ratio of 3:2 are seen for three other dynamically-confirmed black holes (XTE J1550-564, 
GRO J1655-40 and GRS 1915+105). 

About one year after the onset of the 2003 out burst, bipolar X-ra y jets were discovered 
and observed a total of three times using Chandra (ICorbel et al.ll2005l ). Radio observations, 
which commenced several months before the X-ray observations, resulted in four detections 
of t he eastern jet (only ), followed about two months later by a single detection of the western 
jet (ICorbel et al.ll2005l ). Large-scale X-ray jets are rare, having been previously observed for 
only one other microquasar, namely XTE J1550-564 , which is similar in many respects to 
H1743 (for comparisons, see iMcChntock et al.l W0% . In its 1998 outburst, XTE J1550- 
564 produced relativistic jets at early times whose launch date was unambiguously tied to 
the o ccurrence of a remarkable X-ray flare (IHannikainen et al.ll2009l : ISteiner fc McClintock 
2012h . 



H1743's X-fiare on MJD 52766 showed striking similarities to the giant flare of XTE 
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J1550-564 (see Figure 12 of lMcClintock et al.ll2009l ). Of particular note, both flares occurred 
during a dip in the X-ray rms power (0.1-10 Hz), a jump in frequency of the low-frequency 



QPOs, onset of the high-freq uency QPOs, and the apex of power-law emission ( ISobczak et al. 



2000l : iRemillard et al.ll2006[ ). The similar character of these two flares, and the coincidence 
for XTE J1550-564 between the X-ray flare and the launch date of the jets, motivated us to 
search the VLA archive for additional observations of H1743. This search was fruitful, and 
in Section [2] we report three additional radio jet detections at early times that link the jets 
to the 2003 X-ray flare. 

To deduce the source distance and jet inclination angle of H1743, we model the proper- 
motion data derived from the X-ray and radio observations. In doing s o, we closely follow our 



recen t study of the large-scale X-ray/radio je ts of XTE J1550-- 564 (ISteiner &: McC 



intock 



20121 ) , which builds on the pioneering work of IWang et al.l ( 120031 ) and iHao &: Zhang ( 120091 ) . 
Using a model originally applied to gamma-ray bursts, we concluded that XTE J1550-564 
is embedded in a pc-scale cavity in which the jets expanded unimpeded until they impacted 
the cavity walls and rapidly decelerated. We apply this same model to H1743 and obtain 
constraints o n the distance and jet inclin ation angle (presumed to be the inclination of the 



spm axis; see 



Steiner fc McClintockl 120121 ). 



The evolution of H1743's jets have already been studied by lHao &: Zhangl ( 120091 ): how- 
ever, our aims differ from theirs. They were primarily interested in the environment of the 
black hole. While assuming an earlier launch date for the jets, they adopted the nomina l 
values of distance and inchnation {D = 8 kpc and i = 73°) suggested by lCorbel et al.l ( 120051 ). 
Our attention is focused on deriving accurate constraints on D and i for H1743, which we 
use in turn to constrain the spin of the black hol^. 



We measure the spin of H1743 using the continuum-fitting (CP) method ( iMcClintock et al 



We measure the spin 


2011 




Zhane et al. 


1997) 



tion disk Ri^, which is identified with the radius of the innermost stable circular orbit -Risco- 
Knowing both i?isco and M is equivalent to knowing the spin parameter a* because Risco/M 
is a monotonic function o f a^, decreasing from 6 to 1 as the spin parameter increases from 
to 1 dSardeen et "aDll972[ lR. In the CP method, one determines -Risco by modeling the X-ray 



continuum spectrum of the dominant thermal component using a fully relativistic model of 
a thin accretion disk. The observables are X-ray flux, temperature, distance D, inclination 
i, and mass M. In order to obtain reliable values of a^,, it is essential to select X-ray spectra 



^We express black hole spin in the customary way as the dimensionless quantity a* = cJ/GM"^ with 
a* I where M and J are respectively the black hole mass and angular momentum. 

^Using c — G = l. 
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that have a strong thermal component and to have accurate estimates of D, i, and M. For 
H1743, we use our jet model to determine the first two parameters, and we constrain M 
using the known distribution of black hole masses for X-ray transient sources. 



Data 



To search for the presence of radio jets near the time of their expected production 
(Section [1]), we examined high sp atial resolution A-config uration VLA images taken early 
during H1743's 2003 outburst fsee iMcClintock et al.l[2009[ ). Calibrated data from the VLA 
archive for program AR523 on MJD 52779.4, 52782.4, and 52786.4 were imaged using the 
Astronomical Image Processing System (AIPS) task IMAGR. The source was detected at 
8.4 and 14.9 GHz, but here we only use the 14.9 GHz data, which had sufficient angular 
resolution to clearly resolve source components. The synthesized beam was approximately 
0.6" by 0. 2" elongated north -south. Fortunately, the jet position angle is almost exactly 
east-west (ICorbel et al.ll2005l ). allowing us to identify components separated by > 0.2". At 
all three epochs, the source displayed a dominant component and a weak component offset 
towards the west. At MJD 52779.4, just 13 days after H1743's X-ray flare, their separation 
was 166 ± 20 mas. Later, on MJD 52782.4 and MJD 52786.4 the separations were 256 ± 20 
mas and 288 ± 20 mas, respectively. 

The majori t y of t he jet data considered in our analysis are taken from Tables 1 and 3 
of ICorbel et al.l (120051 ). These tables provide jet-source separation measurements for radio 
and X-ray observations which were conducted from 6 months onward following H1743's jet- 
launching flare. The X-ray data consist of three ^ 30 ks Ch andra X-ray observations in 
which both jets were detected. In radio, ICorbel et al.l (120051 ) report on flve observations 
from the Australian Telescope Compact Array (ATCA). The eastern jet was present in each 
image, but the western jet was detected only in the flnal observation. These X-ray and radio 
observations were carried out between MJD 52955 and MJD 53092, when the jet-source 
separations were in the range ~ 4" — 7". The substantially larger angular separations of the 
eastern jet indicate that it is approaching and the western jet is receding. 

In determining the spin of H1743, we analyze the full set of RXTE PCU-2 "standard 2" 
data obtained during the 2003 outburst, with th e spectra binned into 170 half- day intervals. 
These spectra have been modeled in detail by iMcClintock et al.l (120091 ) and ISteiner et al. 
(120091 ) ■ and we use the same data reduction procedures here. Briefly, all the data are 
dead-time corrected, background subtrac ted, and analyzed with the inclusion of a 1% sys- 
tematic uncerta i nty (Jahodaet al.l 120061 ). We standardize all detector calibrations to the 
Toor &: SewardI (Il974j ) values for the Crab using a custom model which adjusts both the 
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overall flux normalization and the spectral shape (see lSteiner et al.ll2010l ). During the early 
weeks of the outburst cycle, RXTE^s pointing was offset by 0.32° from H1743. We have 
corrected the fluxes to the full collimator transmission by assuming a triangular response 
with FWHM = 1° fsee lsteiner et aliboooh . 



The Ballistic Jets: Model and Results 



Our jet model, which is based on one developed by lWang et al.l ((20031), was first applied 
in describing gamma-ray-bursts. Here, we consider a pair of symmetric jets, each ejected 
with an initial kinetic energy Eq and Lorentz factor Tq. During their expansion, the jets 
decelerate as they sweep up gas in their paths. Assuming adiabatic expansion, the evolution 
of each jet is governed by: 



En 



(r - l)Moc2 + a{Tl^ - l)mswc' 



where F is the bulk Lorentz factor of the jet, Mq the mass of the ejecta, a is a numerical 
factor of order unitjifl, and Tsh is the Lorentz factor of randomly accelerated particles at the 
shock front. The entrained mass, nig^, is given by ms„ = Q^nipmiR^ /3, where 9 is the jet 
half opening angle, n the gas density, rrip is the mass of a proton, and R the distance traveled 
by the jet. 

We evolve Eqn. [1] in 2-hour time steps, using the inclination of the jet axis to the 
observer's hne of sight (9) to calculate the projected separation (5) between each jet and the 
central source: 6{t') = R{t)sm 6/D. Here, t' = t ± R{t)cos 6/c is the observer's time, which 
takes into account for each jet the time delay between H1743's rest frame and the frame of 
the observer. 

Our full model requires just five parameters: D, 9, Fq, the launch date Tq, and E, the 
effective energj]^. Because of the association with the X-ray flare, the prior on the launch 
date is taken to be MJD 52766 ± 5 days; we adopt a fiat prior on 9, D, log(Fo), and log{E). 
Our model is fitted via a Markov chain Monte C arlo (MCMC) routine developed using 
the Metropolis-H astings algorithm ( Hastings! Il970[) which has been previously applied with 
this jet model in Isteiner fc McClintockl J20m " The chains are evolved until they are well 



•^cr ranges from 0.35 in the ultrarelativistic lim i t to .73 in the nonrelativistic limit. For additional details 
concerning our model, see lSteiner fc McClintockl (|2012l ). 

I^E = £:o("-/10"^cm-3)-i(e/l°)-2. Following our approach for XTE J1550-564 |Steiner fc McChntock 



2012I ). we scale Q and n using typical values, with density 100 times lower than for the interstellar medium. 
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Fig. 1. — Our best fit model for tlie motion of H1743's radio and X-ray jets. The eastern jet 
is marked by filled circles and the western jet by open circles. Fit residuals are shown in the 
bottom panel using a shght offset in time between eastern and western jets. 
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Fig. 2. — Marginalized probability densities from the MCMC model are shown to arbitrary 
scale. The prior for each parameter is indicated by a dashed line. Fq is constrained by its 
prior at large values, but the other parameters show little dependence on their priors. 
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converged, using ~ 2 x 10^ elements totajj. 

From the VLA data alone, the identification of the pair of radio sources is ambiguous. 
We have applied our model by attributing to the two radio sources each allowed combination 
of eastern jet, western jet, and core. The most probable interpretation is that the two sources 
correspond to emission from the core and the western jet. Alternative pairings are ruled out 
at > 97% confidence by our model. 

The best fit achieved by the MCMC run is shown in Figure [T] and reaches a goodness 
of fit yi^ jv = 4.9/9 = 0.54. Obviously, further modification to the model is not needecj^. 
Distributions for the model parameters are shown in Figure [2l Of chief importance, we find 
that distance and inclination are well constrained: D = 8.5 ± 0.8 kpc and i = 75° ± 3°. 
This distance places H1743 near the Galactic center, which is expected, given its projected 
separation of only ^ 2° from the Galactic center. The time at which the jets were produced is 
constrained to Tq =MJD 52767.6±1.1 days, independent of the prior. This timing supports a 
connection between H1743's X-ray flare and the production of its jets. The speed of the jets, 
Fq, has a relatively low maximum a posteriori estimate, Fq ~ 1.4, but is poorly constrained 
at high values and tracks its prior. For the kinematic energy of each jet, we obtain a large 
uncertainty of ~ 0.5 dex centered around E ^ 10^^ erg. This implies that H1743's jets are 
only about a tenth as energetic as those produced in the 1998 outburst of XTE J1550-564 
or, alternatively, for H1743 either (1) the density of the surrounding medium is much lower 
or (2) the jet opening angle is substantially smaller. 



4. X-ray Continuum-Fitting Analysis 



We now estimate the spin of H1743 by fitting its X-ray spectra. For the three cru- 
cial input parameters, we use the values of D and i derived in the preceding section and 
the distribution of black hole masses discussed bel ow. All of our a nalysis is performed 



using XSPEC vl2.7.0 (lArnaudl Il996l ). Following ISteiner et al.l (120091 ) and making minor 
adjustments, our spectral model has the form tbabs(simpL(8>kerrbb2), where tbabs and 
KERRBB2 are, respectively, the low-energy-absorption and accretion-disk components. The 
component SIMPL scatters a fraction of the thermal disk photons into a Compton power law. 
For H1743, this simple convolution model describes only the broad continuum components 



^An additional ~ lO'^ chain elements, which were generated during training and burn-in phases, were not 
used. 



^We have explored the asymmetric models employed for XTE J1550-564 (jSteiner fc McClintoclal2012r) : i 
and D are unchanged. 
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and is unaffected by the inclusion of weaker features due, e.g., to warm absorbers or spectral 
reflection. 

The four free parameters of the spectral mode0 are the (1) fraction of thermal photons 
/sc scattered into the Compton power law; (2) power-law index F; (3) spin parameter a*; 
and mass accretion rate M. Mass, inclination, and distance are varied in 5 x 10^ Monte-Carlo 
samples, and all 170 spectra are fitted for each setting. Uncertainty in the absolute calibra- 
tion of the X-ray flux is accounted for by randor nly varying the over all flux normalization by 
10% for each triplet setting of M, i, and D (e.g.. ISteiner et al.ll201ll ). We similarly marginal- 
ize over uncertai nty in the viscos i ty parameter a by randomly assigning either a = 0.01 



or a = 0.1 (e.g.. King et al.l 120071 : iPessah et al.l 120071 ). which are two representative values 



available to our model. Both of these uncertainties have a small effect on a^, compared to 
our dominant uncertainty, the unknown black hole mass. (Uncertainty in the mass accounts 
for 50% of our final uncertainty in a^,.) 

For each of the 5 x 10'^ parameter settings, we apply our standard data selection criteria: 
disk luminosity between 3% and 30% of the Eddington lim it: goodness of fit x^/^ < 2; and a 
power-law normalization fsc < 25% (ISteiner et al.ll2009l ). Typically, about 30 spectra pass 
this screening. Finally, each of the 5 x 10^ samples is given a weight according to the mass 
distribution assumed, and random draws are made from the selected spectra to achieve an 
estimate of spin. The dependence between the inferred value of spin and the black hole's 
mass is illustrated in Figure [31 As mass is varied from M = 5Mq to 15M0, spin changes 



from a* 



-0.25 to a, ^ 0.75. 



Recently, lOzel et al.l (120101 ) compiled all the dynamical measurements of mass for black 
hole transients and determined the following best-fit probability distribution, which we 
adopt: 

Exp[(6.30Mq - M)/1.57M0]/1.57M0, M > 6.3Mq, 
0, M < 6.3M0. 



P{M) 



(2) 



In Figure m we show the spin which results when the lOzel et al.l (l2010l ) mass distribution 
is assumed. We find a* = 0.201q;33 (68% confidence interval) with a 90% confidence interval 
of —0.33 < a* < 0.70. We also find that extreme values of spin are ruled out; a^, < 0.92 at 
99.7% confidence. This makes H1743 one of a growing po pulation of black hole microquasars 
known to have moderate spin (e.g., A0620-00, a* 



0.1; iGou et all 12011 XTE J1550-564, 



^Column density A^h is frozen at 2.0x10^^ cm"^ (jBlum et al.ll2009l ). For kerrbb2, limb darkening and 
returning radiation are switched on and the torque at the inner boundary is set to zero. For SIMPL, we use 
the faster upscattering-only option. 
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Fig. 3. — The dependence of spin on black hole mass. These estimates incorporate all sources 
of measurement error. The solid line tracks the average spin at each mass, and the associated 
68% confidence interval corresponds to the region between dashed lines. 

R,,ro IM:2 3 4 56789 
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Fig. 4. — Spin probability for the dimensionless measurement variable RiscojM (top axis) 
which uniquely determines the spin parameter a.^ (bottom axis). This spin estimate is 
obtained by using the adopted transient black hole mass distribution from lOzel et al.l ( 2010 ) 
(solid curve). The vertical solid line indicates the maximum likelihood spin, while the 90% 
confidence range is bounded by dashed verti cal lines. For coni parison, we also show the spin 
estimate from using the mass distribution of iFarr et al.l ( 120111 ) (dotted curve). These results 
take into account uncertainties in M, i, D, a and the absolute X-ray flux calibration. For 
illustration, the distributions over Ri^co/M have been smoothed using a Gaussian kernel 
with 10% width. 
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CS; ISteiner et al.lboill ). 



In Figure HI we also show result s usin g the mass distribution favored by iFarr et al. 



teOllI ). These authors and lOzel et al.l (|2010[ ) used the same black-hole mass data, but Farr 
et al. found that a power-law distribution gave the best fit, with form P{M) oc M~^'^ over 
the mass range 6.IM0 < M < 23M0 (and P{M) = elsewhere). Comparing the Farr et al. 
distribution with our adopted result, we find that the difference is minor: Aa* ~ 0.05. 



5. Conclusions 



We have modeled the proper motion of the radio and X-ray jets of H1743 that were 
launched during an X-ray flare. Based on our purely kinematic model, we obtain firm 
estimates of the source distance, 8.5 ± 0.8 kpc, and the jet inclination angle, 75° ± 3°. 
Using these constraints on D and i, we fitted all 170 X-ray spectra collected during the 2003 
outburst of H1743, applied our data selection criteria, and derived a relationship between spin 
and black hole mass. We then constrained the mass of HI 743 using an analytic distribution 
for transient systems that are similar to H1743, thereby arriving at our final result: a* = 
0.2 ± 0.3 (—0.3 < a* < 0.7 at 90% confidence). Meanwhile, we rule strongly against an 
extreme value of spin: a^, < 0.92 at 99.7% confidence. Two similar microquasars have been 
identified which also produced powerf ul jets while harboring blac k holes with moder ate spins: 
A0620-00 (a* ^ 0.1: [gou et atlboioh and J1550-564 (a* ^ 0.5: Isteiner eTaPboilh . 



This is the first successful application of the X-ray continuum-fitting method that does 
not rely on any dynamical data to place constraints on one or more of the input parameters 
D, M and i - even the orbital period of H1743 is presently unknown! Our constraint on a^, 
can be tightened once a dynamical estimate of mass has been obtained. 
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